Delayed feedback induced directed inertia particle transport in a washboard potential 
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We consider motion of an underdamped Brownian particle in a washboard potential that is subjected 
to an unbiased time-periodic external field. While in the limiting deterministic system in dependence 
of the strength and phase of the external field directed net motion can exist, for a finite temperature 
the net motion averages to zero. Strikingly, with the application of an additional time-delayed 
feedback term directed particle motion can be accomplished persisting up to fairly high levels of 
the thermal noise. In detail, there exist values of the feedback strength and delay time for which 
the feedback term performs oscillations that are phase locked to the time-periodic external field. 
This yields an effective biasing rocking force promoting periods of forward and backward motion 
of distinct duration, and thus directed motion. In terms of phase space dynamics we demonstrate 
that with applied feedback desymmetrization of coexisting attractors takes place leaving the ones 
supporting either positive or negative velocities as the only surviving ones. Moreover, we found 
parameter ranges for which in the presence of thermal noise the directed transport is enhanced 
compared to the noise-less case. 

PACS numbers: 05.40.-a, 02.50.Ey, 05.60.-k, 05.65.+b, 87.15.-v 



I. INTRODUCTION 



Transport phenomena play a fundamental role in many 
physical systems. In this context, the theme of a ratchet 
dynamics has attracted considerable interest over the 
past years. This is particularly due to the fact that the 
ratchet effect assists the creation of a directed flow of par- 
ticles without the presence of any bias force in the sys- 
tem. The ratchet dynamics has been mainly applied to 
biological or mesoscopic systems where Brownian motion 
in a periodic asymmetric potential together with dissipa- 
tion enters in some form the problem and the directed 
motion is generated from nonequilibrium noise, see the 
various overviews on molecular and Brownian motors in 
Refs. [I|-[U. On the other hand for periodic systems 
with maintained spatial symmetry the accomplishment 
of directed net motion necessitates that the system is 
exerted to additional biasing (symmetry-reducing) im- 
pacts. Here we consider particle motion in a washboard 
potential which is often employed paradigm to 

model transport in one-dimensional periodic and sym- 
metric structures @|,0]-O- A symmetric (unbiased) 
external force is assumed to rock the potential. Our 
aim is to demonstrate that directed particle transport 
can be achieved with the application of a time-delayed 
feedback method in a wide temperature range. Although 
the delayed feedback method was originally proposed by 
Pyragas [l4[ to stabilize unstable states in determinis- 
tic systems meanwhile it has been facilitated in various 
other contexts [l5[ among them there is also the control 
of purely noise-induced oscillations [lj| , [I?) ■ Recently in 
the context of controlling transport in Brownian motors 
a feedback strategy has been successfully utilized for two 
ratchet systems interacting through a unidirectional de- 



lay coupling [18| ■ The effect of time-delayed feedback on 
the rectification of thermal motion of Brownian particles 
has been studied in overdamped ratchet systems [19j]-[25|. 
A recent experimental implementation using such a feed- 
back mechanism for a flashing ratchet has been realized 
with an optical line trap: it has been observed that the 
use of feedback increases the ratchet velocity up to an 
order of magnitude [2t| . in agreement with theory. 

Stabilization of chaotic motion in deterministic iner- 
tia ratchet systems to increase the current efficiency was 
considered in Refs. 27|,[28[. Furthermore, an asymmet- 



ric ratchet potential with included time-delayed feedback 
was treated in the context of an inertial Brownian motor 

Our paper, dealing with time-delayed feedback induced 
directed motion, is organized as follows: First we intro- 
duce the model of an inertial Brownian particle evolv- 
ing in a symmetric spatially periodic potential under the 
influence of an additional time-delayed feedback term. 
The subsequent section concerns the underlying deter- 
ministic dynamics. In particular bifurcation diagrams 
with and without applied time-delayed feedback are dis- 
cussed. The impact of a heat bath of fixed temperature 
on the particle transport features is studied in Section 
IIVI Finally we summarize our results. 



II. THE MODEL 



We consider an inertia Brownian particle that is mov- 
ing along a one-dimensional periodic structure. The dy- 
namics is governed by the following inertial Langevin 
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equation expressed in dimensionless form 

q + jq = ~-r + F sin(w t + 6 ) + £(*) + f(t). (1) 
aq 

The dot denotes differentiation with respect to time. The 
particle evolves in a spatially-periodic and symmetric po- 
tential 

U(q)= C/(g + l) = -cos(27rg)/(27r), (2) 

of unit period L = 1 and barrier height AE — 1 jix and its 
position and velocity are quantified by the variable q{t) 
and q{t) = respectively. The particle is driven by 
an external, time-dependent forcing field of amplitude F, 
frequency u> and phase 9q- I n addition it is subjected to a 
Gaussian distributed thermal, white noise £(t) of vanish- 
ing mean (£(£)) = 0, obeying the well-known fluctuation- 
dissipation relation (£(*)£(*')) = 2jk B T5(t - 1') with k B 
and T denoting the Boltzmann constant and tempera- 
ture, respectively. The friction strength is measured by 
the parameter 7. The last term in Eq. ((T|) denotes a 
continuous time-delayed feedback term of the form 

f(t) = K[q(t - t) - q{t)\ (3) 

of strength K and with delay time r. Before we embark 
on the study of the Brownian particle motion it illustra- 
tive to consider the deterministic limiting case arising for 
vanishing thermal noise, i.e. T = 0. For our study we fix 
the following parameter values 7 = 0.1, uj = 2.25, 8 = 0. 
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FIG. 1: (a): Bifurcation diagram as a function of the ampli- 
tude of the external driving without time-delayed feedback, 
i.e. K = and remaining parameter values: u> = 2.25, 7 = 0.1 
and T = 0. (b): Mean velocity v m /vo as a function of the am- 
plitude of the external driving. For later use a dashed vertical 
line is drawn at the value F = 1.5 



III. THE DETERMINISTIC CASE 

The dynamics of the deterministic system (T = 0) ex- 
hibits very rich and complex behavior and depending on 
the parameter values and initial conditions one finds pe- 
riodic, or aperiodic (quasiperiodic and/or chaotic) solu- 
tions in the long time limit [3l|-[3g|. The character of the 
phase flow evolving without feedback term (K = 0) in a 
three-dimensional phase space is conveniently displayed 
by a Poincare map using the period of the external force, 
T e = 2n/u! = 27r/2.25 ~ 2.791, as the stroboscopic time. 
The deterministic equation of motion was integrated nu- 
merically and omitting a transient phase points were set 
in the map at times being multiples of the period dura- 
tion T e . In Fig.[T](a) the bifurcation diagram as a func- 
tion of the amplitude of the external driving is depicted. 

Particle transport is quantitatively assessed by the 
mean velocity which we define as the time average of 
the ensemble averaged velocity, i.e. 

«m = i / " dt'{v n {t')), (4) 
J s JO 

with simulation time T s and with the ensemble average 
given by 

1 N 

(Vn(t)) = W J2 V "( t )- ( 5 ) 



Here N denotes the number of particles constituting 
the ensemble with associated initial conditions gVi(0) and 
iVi(0) that are chosen such that the interior of the unper- 
turbed separatrix in phase plane is uniformly covered. 
We express v m in terms of the ratio of the spatial and 
temporal periods L/T e = vq with vq — 0.358 being the 
velocity for running solutions that advance by one spatial 
period during one period duration of the external field. 

For undercritical amplitudes of the modulation force 
F < 0.6 those particles which are initially residing near 
the bottom of a potential well, remain trapped. Increas- 
ing F leads to escape from the potential wells and the 
particle jumps subsequently from one well to another 
one. The arising two typical scenarios are the pinned 
and running states respectively. In the former state the 
motion proceeds at most over a finite number of spa- 
tial periods whereas in the latter state motion is directed 
and unrestricted in the spatial dimension. In terms of 
the phase flow running asymptotic solutions correspond 
to phase locked periodic attractors transporting a parti- 
cle with velocity v = m/n over m spatial periods of the 
potential during n period durations T e of the external 
periodic field. Running asymptotic solutions may also 
be supported by aperiodic attractors. 

In the bifurcation diagram associated with the dy- 
namics without applied time-delayed feedback shown in 
Fig.[Tj(a) one recognizes vertically extended stripes cov- 
ered densely with points corresponding to non phase 
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FIG. 2: As in Fig. [T] but with switched-on time-delayed feed- 
back of strength K — 0.9 and delay time r = 1.95. The 
remaining parameter values read as: u = 2.25, 7 = 0.1 and 
T = 0. 



locked aperiodic attractors and several periodic windows 
as well as period-doubling cascades to chaos. These fea- 
tures of the phase flow are readily attributed to the 
resulting mean velocity of the net motion (depicted in 
Fig. [U(b)). The ensemble average is taken over an ensem- 
ble of N = 5000 trajectories with uniformly distributed 
initial conditions q(0) and i>(0). For computation of the 
long-time average the simulation time interval for each 
trajectory is taken as T s — 5xl0 5 ~ 1.8xl0 5 xT e . We no- 
tice almost in the entire .F-range vanishing mean velocity 
v m = 0. The exceptions are the intervals 0.74 < F < 0.80 
and 1.36 < F < 1.62 for which the solutions are associ- 
ated with multiple coexisting attractors lying in fairly 
extended periodic windows in Fig. QJa). Focusing inter- 
est on the latter one we note that at F > 1.36 tangent 
bifurcations give birth to two coexisting period-one at- 
tractors. The upper one of them is related with positive 
particle velocity v = vo > whereas on the lower one 
particles move with velocity v = —vq < 0. For increasing 
F > 1.62 these attractors are destroyed by way of crisis 
after passage through a period-doubling route to chaos. 

The oppositely running solutions attributed to the two 
period-one attractors contribute to the mean velocity 
with different weight with the one with positive velocity 
v = vq being dominant and thus yielding the window of 
positive mean velocity v m . Apparently, the directed mo- 
tion results from a lowering of the dynamical symmetry 
caused by the external modulation field [3(3], [l3[. That 
is, even though the potential and the external modula- 
tion field are spatially and temporarily symmetric respec- 



tively, with the choice of a fixed phase #0 the symmetry 
of the flow is reduced and a phase-dependent net motion 
is found. (Note that additional averaging over the phase 
#0 yields zero mean velocity.) Due to symmetry reasons 
it holds that the sign of the mean velocity is reversed 
upon the changes 9q — — > 60 = tt and F — > —F re- 
spectively. However, there exists a phase < 8q < ir for 
which symmetry between the two coexisting periodic at- 
tractors supporting solutions with velocities of opposite 
sign, vo and — vq, is restored and therefore the net motion 
vanishes. 

In the numerical simulation of the system |T]) with ap- 
plied time-delayed feedback term ([3]) we set f(t) =0 in 
the interval t 6 [0,r), that is the system is affected by 
f(t) only for t ^ r. We performed extensive numerical 
studies to identify optimal parameters of the feedback 
term which establish efficient directed net motion. It 
turns out that this is achievable for delay times in the 
range of 1.65 < r < 2.00 and for feedback strength K > 
0.8 (see also further below in Fig. [6]). In the following we 
illustrate exemplarily the impact of time-delayed feed- 
back on the transport properties for a feedback strength 
K = 0.9 and delay time t = 1.95 ~ 0.7 x T e . With 
such appropriate feedback term applied the extension 
of the aperiodic regions shrinks considerably and only 
a comparatively narrow band of aperiodic behavior for 
1 < F < 1.11 prevails in the bifurcation diagram illus- 
trated in Fig. [2](a). Remarkably, for F > 1.63 the lower 
period-one attractor supporting negative velocities looses 
stability and is converted into a repellor leaving the up- 
per attractor of positive velocity as the only persisting 
attractor. Thus application of feedback results in a re- 
shaping of the bifurcation diagram. In fact, due to the 
absence of the period-one attractor supporting motion 
with velocity v = — vq only running solutions with veloc- 
ity v m = vq are then recognized in Fig. [2](b) showing the 
mean velocity as a function of F. Otherwise the mean 
velocity raises form zero level for overcritical F > 1.1 and 
grows with increasing F until F > 1.63 when v — vq is 
attained. 

In order to gain insight into the feedback-induced 
mechanism that leads to directed transport (occurring for 
F > 1.63 in the period-one window in Fig. [2]) we display 
the temporal behavior of the feedback term f(t), given in 
Eq. ([3]), and the external driving term F(t) — Fsm(ojt) 
and their sum F(t) + f(t) in Fig. [3] for driving ampli- 
tude F — 1.8. Throughout the time the feedback term 
f(t) performs oscillations possessing considerable asym- 
metry. Most importantly, these oscillations are entrained 
to the (symmetric) external driving term with a phase 
shift. The sum F(t) + f(t), performing asymmetric os- 
cillations, determines the effective rocking force exerted 
on the particle. This ratching force is self-induced due to 
the feedback in comparison with the externally imposed 
ratching force in the form of asymmetric periodic driving 
fields considered in 30]. 

How the directed rightward particle motion is enforced 
by this effective biasing rocking force is illustrated in 
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FIG. 3: Time evolution of the feedback term f(t), the driving 
force term F(t) and their sum F(t) + f(t) with K = 0.9 and 
delay time r = 1.95. The remaining parameter values read 
as: F = 1.8, w = 2.25, 7 = 0.1 and T = 0. 
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FIG. 4: Time evolution of the particle amplitude q(t) and the 
effective rocking force F(t) + f(t) with K — 0.9 and delay time 
r = 1.95. The remaining parameter values read as: F = 1.8, 
lj — 2.25, 7 = 0.1 and T = 0. Horizontal solid (dashed) lines 
indicate the position of the maxima (minima) of the unbiased 
washboard potential. 



Fig. [U (For the present discussion we discard the contri- 
bution from —dU(q)/dq to the total force.) Clearly, an 
oscillating rocking force leads to passages of forward and 
backward motion, also called enhancement and depreci- 
ation periods. Crucially, the feedback term is suitably 
entrained to the external driving term in such a way that 
the period of forward motion is longer than its backward 
counterpart, denoted by Tf and Tj, respectively in Fig. 2] 
To be precise, at the moment when the backward motion 
of the particle terminates the effective rocking force has 
not yet reached its maximum as indicated by the left ver- 
tical line in Fig. 3J Subsequently the particle moves in 
the forward direction due to the ongoing positive rock- 



ing force that after passing through its maximal value 
declines. Nevertheless, until half of the time span Tf is 
reached the particle motion is still enhanced in the for- 
ward direction. Afterwards for t > Tj/2 the effective 
rocking force becomes negative and thus the momentum 
of the particle is reduced steadily with increasingly nega- 
tive values of F(t) + f(t). The end of the forward-motion 
period is designated by the middle vertical line in Fig. 2] 
However, there remains only comparatively little time, 
namely Tb < Tf, during which backward-motion is en- 
forced, that is when during the depreciation period the 
effective rocking force is negative. 

Consequently, the asymmetry in the enhancement 
and depreciation phases serves for a rather long period 
of forward-motion compared to the backward-motion. 
Therefore the effective motion of the particle proceeds 
to the right. Notably, this feature is induced by entrain- 
ment of the asymmetric time-delayed feedback term to 
the symmetric external modulation field if the feedback 
strength and delay time are suitably chosen. Notice that 
this oscillation behavior of the time-delayed feedback 
term is different from noninvasive control methods where 
the delayed feedback control vanishes once a targeted un- 
stable periodic orbit has been stabilized (l4T| . flol] . [28| . [29| . 



IV. DIRECTED THERMAL NET PARTICLE 
MOTION 

We now study the impact of a heat bath of fixed tem- 
perature T > on the particle transport features. With 
the inclusion of finite thermal noise transitions between 
the now metastable attractors are likely and, indepen- 
dent of the initial conditions, trajectories permeate the 
whole phase space rendering the dynamics ergodic. 

For the computation of the mean velocity the ensemble 
average in © was taken over N = 5000 realizations of 
the thermal noise for an arbitrarily chosen initial condi- 
tion. The mean velocity v m represented without feedback 
term but in the presence of thermal noise of a fairly high 
level of ksT = 0.1 x AE in Fig. E](a) is zero regard- 
less of the value of the modulation field strength F. This 
has to be distinguished from the preceding noise-less case 
where there exist even for K = regions of nonvanishing 
net motion (cf. Fig. [I} due to the fact that motion on 
attractors with different sign of the velocity contribute 
with distinct weight to the asymptotic net motion. In 
other words, the impact of the noise leads to symmetriza- 
tion of the basins of attraction of transporting periodic 
and/or aperiodic attractors. Therefore any initial condi- 
tion yields zero asymptotic current. 

Interestingly, this situation changes imposing the 
Langevin dynamics additionally to the time-delayed feed- 
back and we found parameter constellations for which di- 
rected net motion results despite the presence of strong 
noise. 

As Fig. H](a) reveals, applying feedback of strength 
K = 0.9 and delay time r = 1.95, the mean velocity 
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FIG. 5: (a): Mean velocity in dependence of the external 
modulation field amplitude without and with applied time- 
delayed feedback force as indicated in the plot and for thermal 
energy ksT — 0.1 x AE. (b): Mean velocity in dependence 
of the thermal energy expressed in units of the barrier energy, 
ksT/AE, and for fixed F = 1.5 and delay time r = 1.95. 
The remaining parameter values are lu — 2.25 and 7 = 0.1. 



as a function of the external modulation field strengths 
exhibits a resonance-like structure for 1.2 < F < 1.7, 
i.e. for values of F for which transport exists in the de- 
terministic case. Strikingly, with the impact of thermal 
noise the feedback-controlled transport proceeds more ef- 
ficient in comparison with the deterministic case in the 
fairly wide range 1.2 < F < 1.7 (cf. Figs. M and [5](a)). 
On the other hand, for F > 1.7 the directed transport 
feature of the deterministic system is destroyed by the 
thermal fluctuations. Furthermore, like in the determin- 
istic case, there exist a threshold value for the feedback 
strength beyond which directed net motion is achieved 
(see Fig. [6](a)) and the range of delay times being opti- 
mal for running solutions is indicated by the resonance- 
like structure in Fig.[5](b). 

^From Fig. [5](b), illustrating the mean velocity in de- 
pendence of the noise strength, one infers that directed 
transport is sustained up to comparatively high noise 
strength before it ceases eventually to exist. On the other 
hand, it is seen that without feedback term, i.e. K = 0, 
the mean velocity averages to zero already for low noise 
intensity. As the damping strength is concerned we found 
that feed-back-induced directed motion is maintained as 
long as the system remains in the underdamped regime, 
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FIG. 6: (a) Mean velocity in dependence of the strength of 
the time-delayed feedback and for fixed F — 1.5 and r = 1.95. 
(b): Mean velocity in dependence of the delay time and for 
fixed F = 1.5 and K — 0.9. The remaining parameter values 
are u = 2.25, 7 = 0.1 and k B T = 0.1 x AE. 

that is for 7 < 2.5 ~ ljo with luq = \^2tt being the fre- 
quency of harmonic oscillations near the bottom of a po- 
tential well. 

Further insight into the origin of running solutions 
in the noise case is gained from Fig. [7] showing the 
stroboscopic map of the average trajectory defined as 

la = S„=i Qn/N and v a = X) n =l Vn l N with and witl> 
out presence of the time-delayed feedback for external 
modulation strength F = 1.5, a value for which in the 
deterministic case the mean velocity remains nearly un- 
affected when the feedback is applied (to ease the eyes 
the vertical dashed line is drawn at the position of the 
value F = 1.5 in Figs.[j]2]and[5](a)). Without feedback, 
K = 0, the velocity range of the average trajectory is not 
only symmetric with respect to v a — but remains con- 
fined within the boundaries of the separatrix loop of the 
undrivcn, undamped deterministic system (7 = F = 0). 
Thus the solutions represent pinned states. With applied 
feedback of strength K = 0.9 the v-symmetry is broken. 
Moreover, the stroboscopic plot of the average trajectory 
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FIG. 7: Stroboscopic map of the average trajectory without 
(K — 0) and with applied feedback (K = 0.9 and r = 1.95) as 
indicated in the plot. The remaining parameter values read 
as F = 1.5, r = 1.95 and k B T = 0.1 x AE. The attractor 
with v = vo > of the deterministic system arising with the 
applied feedback is also shown. The dashed lines indicate the 
separatrix of the deterministic conservative undriven system, 
i.e. 7 = F = 0. We also superimposed the projection of the 
orbit of the deterministic system on the v — q— plane (blue 
line) . 



densely covers the curve obtained when the periodic os- 
cillations of the noise-less system are projected onto the 
q — v— plane. This implies the existence of an attractive 
curve related to near-torus-motion in phase space for the 
noise case. The corresponding stroboscopic map of the 
dominant deterministic periodic attractor with velocity 
v = vq is also drawn in Fig. [7J Crucially for directed net 
motion with v m = vq arising for T > the average tra- 
jectory sticks to the near-torus motion never exploring 
other parts of the phase space during the whole simu- 
lation time interval T s = 10 5 . The time evolution of 
the corresponding average coordinate q a is depicted in 
Fig. [H Conclusively, despite the fact that the resulting 
motion in the presence of noise and feedback evolves no 
longer perfectly synchronous but still phase locked with 
the external periodic modulation it nevertheless exhibits 
behavior reminiscent of that found in the corresponding 
limiting deterministic system. To be precise the peri- 
odic oscillations of the deterministic dynamics are re- 
placed by near-torus motion in the Langcvin dynamics 
accomplishing phase locked aperiodic transport in the 
sense that on average during one period duration of the 
external field particles move by one spatial period, i.e. 
v m = v = L/T e . 
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FIG. 8: Time evolution of the average coordinate q a (wiggled 
line). The dashed straight possesses slope Do = L/T e ~ 0.358. 
The bullets represent the position of the particle at moments 
being multiples of the period duration of the external modula- 
tion T e — 2-k/u. The inset displays the long-time evolution of 
q a . The parameter values read as F — 1.5, K = 0.9. r = 1.95 
and k B T = 0.1 x AE. 



V. SUMMARY 

In conclusion, we have identified and characterized a 
transport regime for underdamped Brownian particles 
evolving in a symmetric washboard potential under the 
mutual impact of an unbiased external periodic field 
and time-delayed feedback. We have studied first the 
deterministic case of zero temperature. It is has been 
shown that without feedback in some ranges of the am- 
plitude and certain phases of the external modulation 
field windows of directed particle current exist. This 
is due to the fact that attractors associated with phase 
locked oppositely-running solutions contribute with dif- 
ferent weight to the net current. The direction of the 
net current can be reversed with a suitable choice of the 
phase of the external modulation field. 

Interestingly with an applied time-delayed feedback in 
the deterministic system there exist parameter ranges for 
which the attractors supporting velocities of a definite 
sign loose stability and are converted into repellors leav- 
ing the attractors of opposite velocity as the only persist- 
ing ones. Hence, the efficiency of the net particle current 
is improved. There exist values of the feedback strength 
and delay time for which the feedback term performs os- 
cillations that are phase locked to the time-periodic ex- 
ternal field. This yields an effective biasing rocking force, 
promoting periods of forward and backward motion of 
distinct duration, and thus directed motion. 

On the other hand, for a finite temperature and with- 
out time-delayed feedback the net motion averages to 
zero. That is, the symmetry between the coexisting at- 
tractors supporting negative and positive velocities is re- 
stored by the impacting thermal fluctuations. Strikingly, 
in contrast to the case without feedback, with applied 
time-delayed feedback of appropriately chosen strength 
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and delay time we find in a wide temperature range 
complete desymmetrization of coexisting attractors sup- 
porting oppositely running solutions. As a consequence 
the particle motion proceeds exclusively in one direction. 
Therefore the feedback-induced transport is not only ro- 
bust with respect to thermal noise but moreover, there 
exist parameter ranges for which in the presence of ther- 
mal noise the transport is more efficient than in the cor- 
responding deterministic limiting case. We have iden- 
tified attracting curves in phase space which are linked 
with motion on a torus with small deviations supporting 
transport. 

We stress the difference between the transport control 



mechanism described above and the control of transport 
properties in inertia ratchet systems [13] , [HI , [HI facili- 
tating the stabilization of certain targeted (unstable) pe- 
riodic orbits via time-delayed feedback as well as the fea- 
tures in feedback flashing ratchets when the potential is 
alternatively switched on and off in dependence of the 
state of the system [h|-[24]. In our case control of trans- 
port is achieved if the time-delayed feedback term with 
appropriately chosen strength and delay time entrains to 
the external time-periodic field yielding effectively a bi- 
asing force that rocks the washboard potential such that 
distinct durations of the periods of forward and backward 
motion ensue. 
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